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Abstract 

The method of the real time perturbative calculations of nonequilibrium averages is 
generalised to the case of varying chemical potential. Calculations are performed in the 
frame of Zubarev's nonequilibrium density matrix approach. In this approach perturbations 
of temperature and other thermodynamical parameters are taken into account explicitly 
including nonlinear terms. It differs from the Schwinger-Keldysh approach through the 
choice of more general initial conditions for the density matrix. 
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1 Introduction 



Quantum field theory (QFT) at finite temperature in the equilibrium |l|, |], || f|, ^ and 
nonequilibrium || [7], ||, |S| cases has many applications in study of the early Universe 



|TT| , |T2|| and heavy ion collisions |13[. Both kind of systems are intrinsically nonequilibrium 
in nature, and only some details of their behaviour may be considered in the frame of 
equilibrium description. 

The purpose of this paper is to introduce a Green function generating functional for 
the nonequilibrium QFT at finite temperature and chemical potential in the case when 
these thermodynamical parameters depend on space and time coordinates. 

The approach is based on Zubarev's nonequilibrium density matrix []T4] (called by 
author nonequilibrium statistical operator (NSO) [15|]). The NSO is a solution of the 
Liouville equation with initial conditions in the form of a local-equilibrium density matrix. 
This matrix does not describe irreversible processes such as thermal conductivity. It 
corresponds to the maximum entropy of a system for a given distribution of energy, 
charge density etc., but it is not a solution of the Liouville equation, that rules the system 
dynamics. This deficiency plays no role if the initial condition in the form of a local- 
equilibrium density matrix is fixed at the time t — > — oo. Hydrodynamical equations 
which account for the dissipative terms are derived from the NSO explicitly (see [15||). 



A situation with varying temperature is already analized in [HJ and account of a 
chemical potential will be performed here in the same manner. However, there are some 
technical differences. It is shown in |TB[ that consideration of varying temperature requires 
the introduction of a more complicate time integration contour in the complex time in- 
tegration plane in comparison with the equilibrium case @]. It will be shown here that 
in the case of a uniform temperature distribution and varying chemical potential a com- 
plicated contour is not required. Thus, temperature plays a special role in the quantum 
field theoretical description. 

There are other approaches to the study of QFT in inhomogeneous media. It was 
argued in |T7|], Zubarev's f|TJ], [T5j and McLeennan's [JTS], [HJ approaches are the most 
useful, since they take into account thermal and charge inhomogeneities explicitly. These 
approaches have the same numerical results. 

The most used method in QFT, that due to Schwinger-Keldysh 0] is also based on 
the Liouville equation but with equilibrium initial condition for the density matrix. Only 
those thermal perturbations that occur after strong external action on the initially equi- 
librium system are taken into account by this approach . It means that an inhomogeneous 
distribution of the initial values of thermal parameters is not taken into consideration. 
And in general, this approach does not determine an explicit expression for the tempera- 
ture chemical potential and other thermodynamical parameters. 



2 Density matrix 

It is possible to consider nonequilibrium stationary systems in the frame of a generalised 
Gibbs method. The idea of local equilibrium is very useful for this aim. It is based on 
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the relaxation time property of nonequilibrium systems. With this property, the system 
rapidly achieves a local equilibrium state and then slowly relaxes to equilibrium [ f2~T|j . 
Investigation of the second long stage is possible in terms of NSO. As it is shown in [ |i"5|| , 
the following local equilibrium density matrix corresponds to the maximum of the system 
entropy 

P t (4\ _ C-l e -/ d 3 x^"(x,t)f 0l ,(x,t)+/ d 3 xa(x,i)J (x,t) ^ 1) 

where 

(3 v (x) = f3(x)u u (x), a(x) = p(x)ii{x), (2.2) 

(3{x) = 1/T(x) - - inverse local temperature, u u (x) - - hydrodynamical velocity, p{x) 
- chemical potential corresponding to the charge density Jo (p k and Jq if there are 
some conserved charges in the system), T^ v is the energy momentum tensor and S/ - 
the partition function that normalises ( |2.1|) . It is clear that in the limit of uniform 
values of temperature and chemical potential distribution we recover the well-known grand 
canonical Gibbs density matrix, that in the rest frame (where = (1, 0, 0, 0)) is 

Po = H^e-^+^o, (2.3) 
where H is the Hamiltonian of the system and 

Qo= / rf 3 xJ (x,t) (2.4) 



its charge. 

The density matrix ( |2.1| ) gives the correct results for the energy and charge density. 
But it doesn't describe dissipative processes such as thermal conductivity, diffusion, vis- 
cosity etc. 

A density matrix that completely describes the dynamics of the system including dis- 
sipative processes is obtained by Zubarev as a retarded solution of the Liouville equation 



m 



^-i\p,A\ = (2.5) 



with the initial condition in the form of a local-equilibrium density matrix (|2.1|) . For this 
aim an external source that breaks time invariance of the equation is introduced in the 
right hand side of (|2.5|) 

^- l [p £ ,H] = -e(p £ -p l ). (2.6) 

The limit e —>■ is to be taken after thermodynamical averaging. The same procedure of 
the retarded solution selection of the equation of motion with the free initial conditions 
is used in scattering theory ||22|| . The solution of (|2.6|) is [14], [15 



p £ (t)=e dhe^piit + h). (2.7) 
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Thus, integration of (5T6) with the positive sign of e gives the retarded solution of the 
Liouville equation Q2.5Q . This solution corresponds to the entropy 

S = -Tr{p £ \np £ } (2.8) 

increase. There are some other useful techniques to obtain the smooth density matrix 
describing entropy increase (see for example |T5|). The exact solutions p exac t of ( |2~5| ) 
corresponds to the conserved entropy S = — Ti{p exact \np exact }. It may be averaged in 
quantum states with close energies or in time intervals to obtain coarse grained 

Pr(t) = -^^ Ar p exact (t) (2.9) 

or time averaged 

p T (t) = ^ dt' p exact {t') (2.10) 

density matrices. So, St = — Tr{prlnpr} and S T = — Tr{p T lnp T } increase in time. 
In the general case, matrix ( |2.1| ) is expressed 

Pi 



s -l e -/dcrf 1 (/3^(x,t)f M 4x,t)-a(x,t)J M (x,t)) ^ 



where defines some space-like hypersurface. Now it is possible to introduce transformed 
energy-momentum tensor depending on chemical potential. 

n — f da' 1 (/3^(x,t)f MI ,(x,t)-/ i (x,t)/3 , '(x,t)u l ,(x,t)J M (x,t)) _ 

Pi — "Z e J — 

= E^e'I d*»{r(x,t)f^( x ,t)^ (2.12) 

where 

T'^ix) = ffwix) - p{x)u v (x)J^(x), (2.13) 

and identification 

/3(x)=p v (x)u v (x) (2.14) 

was used. 

When the distance of the hydrodynamical parameter variations L is bigger than the 
correlation lengths 

L>/ corr (2.15) 



the hypersurface may be chosen in an arbitrary way |23|] and it is possible to consider 
the hypersurface for which dx° = (da^ 1 = (d 3 x, 0, 0, 0)). In this case NSO may be written 
in the form (127 



Pe 



e f rftlSr i e /rf 3 x^(x, i+il )f^(x, i+il)) (216) 

J — oo 



with Tq v defined by (|2~T3]). 
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3 Thermal perturbations 



Consider the inverse temperature distribution slightly differing from a certain constant 
value Pq = PqUq by a small parameter /3f (x) 



(3.1; 



Using this expression it is possible to apply our results [16], [L7| but with the 
transformed energy- momentum tensor ( |2.13| ) in the NSO ( 2.16 ). NSO may be written in 
the rest frame (uq = (1, 0, 0, 0)). In this frame PqP^ = /3 H 



p £ {t) = e dtie 



_eti"- 1 



cxp 



-PoH'- 



- I d- 3 x/3f(x,t + t 1 )T^(x ) t + t 1 ) 



where 



H' = H 



d 3 xfi(x)J°(x). 



(3.2) 



(3.3) 



Now we make use of a well-known formula for the exponent of noncommuting operator 
sums 



26f rewritten for our case in [[25 
exp \-PqH' + B(r 



c i exp 



c 



dx /3 1 B(t + x ) 



(3.4) 



where C is an arbitrary contour in the complex time-plane appropriate for the real time 
QFT U, |, |7|, ||, H and T c — is an ordering by the time variable along C. With the 
help of (|3.4j) , NSO may be presented as 



xT c < exp 



c 



p e (t)=e I dhe^Z^e-W x 
f^(x,a;o + t + ti) 



/3o 



(3.5) 



Here the chemical potential is to be considered. After a reduction ( |3.3|) , term with the 
chemical potential became part of the Hamiltonian H'. In analogy with ( |3.1|) , it may be 
presented as 

/i(x) = (jl + Hi(x) (3.6) 

and the part of H' dependent on \i\ will be considered as a perturbation. 

Perturbations of thermodynamical parameters are called "thermal perturbations" pT 



differing from "mechanical" ones which may be considered as part of the Hamiltonian. 
However after reduction ( |3.3| ) of the Hamiltonian, perturbations of the chemical potential 
may be formally taken into account as mechanical ones. Analogous technical tricks are 
not possible with the temperature and it leads to the more complicated contour for the 
time integration [p], |TT|| . 
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Thus, instead of an ordinary perturbation in the Hamiltonian V (usually characterised 
by the coupling constant) we have a more complicated one 



V -> V - J d 3 xpi(x)j°(x) 



(3.7) 



dependent on both coupling constant and chemical potential deviation ni(x). The first 
term in (|3.7|) corresponds to the familiar vertex and the second one will be presented in 
the diagram as an insertion into the line. 



4 Generating functional of the theory 

Now a generating functional for the Green's functions may be considered 



Z\j] = Tr {p e Te>*} = e j" dhe^N^Tr { Po T c 



J3iT> 



where 



J4> 



cxp 



i / d A yj{y)<P{y) 



cxp 



— i I d 4 x 

c 



00 



To^xo + t + h) 



and N(ti) is the residual part from the partition function Sj. 
Using the identification 

r a 

F(4>)e^ = F 



iSj 



i<t> 



we may rewrite (|4.1|) in the form 



x exp 



-ii d x 
lc 



4 /^(x,t+t x ; 



/?0 
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i6j(pc, xo + t + tx) 



/ dt x e etl x 

J— oo 



i'=o 



where T' [S/iSj] is the same function of S/iSj as T'(<f>) of and 

Z[7 ) j1=1t{ATo[^T[^]} 



(4.1) 

(4.2) 
(4.3) 



(4.4) 



(4.5) 



(4.6) 



The generating functional ( |4.6|) may be defined with the unit time-ordering procedure 
if the time-integration contour C (fig. 1) is introduced. This contour differs from the 
one considered in the previous papers ]16|, however it satisfies all the conditions used 
there. Segment C\ is one where external fields fixed in the second exponent in ( |4.6|) and C2" 
and C3 — for thermal perturbations in the first time-ordered exponent in the expression 
( |4.6D . As previously described, the vertical part is considered at —00 and, therefore, it 
plays no role. A difference between the contour indices 2' and 2" will be observable from 
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the region of time integration in the generating functional and it is not specified by the 
separate index value of <p a and J a . This contour consists only of 3 horizontal segments and 



the theory is based on the 3x3 free propagators (instead 4x4 used previously [16. [T 
It is clear that the generating functional has the same form as in ||16|| : 



Z\juh,h] = e f dt ie £h exp{-i f d^xj^{-l) a V 



x exp < i 




dxr, / d x 



o=l 
3 



00 



E(-i)% 



t i 



a=2 



X 



exp<-i / d 4 x^2(-l) a [i 1 (-x, t + ti)J 

{ a=l 
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x + t + £i) 
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i5j a (x, 


x + t + £i) 
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X + t + ti) 



X 



X 



X 



x (4.7) 



6XP I ~\ I d ^ Zld ^ Z2 S 3a{zi)Vf oflo { Zl - z 2 )j b (z 2 ) 



a,b=l 



The difference with the generating functional derived in consists in the presence of the 
chemical potential perturbations fii and in the free Green function that depends not only 
on the constant part of the temperature /3q, but also on the constant part of the chemical 
potential /i Free propagators dependent on the chemical potential are well-known in the 
equilibrium case (see for example |27j ) 

Finally I would like to mention that in the case of a constant temperature and varying 
chemical potential term, (3\ is absent in ( |3.5|) and no double-time ordered product occur 
in ( [4.11 ). It makes it possible to consider a theory in this particular case with the ordinary 
contour that is suitable in the equilibrium case, where free propagators are 2x2 matrices. 

It is important to mention also that this theory is free from pinch singularities. Argu- 
ments for this fact are the same as in |[T6|, IT 
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Figure 1: Contour for the nonequilibrium case 
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